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PACS 6 8 . 3 5 . Rh - Phase transitions and critical phenomena 

Abstract. - We have measured the center-of-mass fluctuations of the height of a contact line at depinning 
for two different systems: liquid hydrogen on a rough cesium substrate and isopropanol on a silicon wafer 
grafted with silanized patches. The contact line is subject to a confining quadratic well, provided by gravity. 
From the second cumulant of the height fluctuations, we measure the renormalized disorder correlator A(u), 
predicted by the Functional RG theory to attain a fixed point, as soon as the capillary length is large compared 
to the Larkin length set by the microscopic disorder. The experiments are consistent with the asymptotic 
form for A (it) predicted by Functional RG, including a linear cusp at u — 0. The observed small deviations 
could be used as a probe of the underlying physical processes. The third moment, as well as avalanche-size 
distributions are measured and compared to predictions from Functional RG. 



o 
o 



> 
in 



o 

ON 

o 

> 

X 



A direct measurement of the fixed-point function A(u), the 
so-called renormalized disorder correlator, which plays a cen- 
tral role in the Functional RG theory (FRG) of pinned elastic 
systems, was recently proposed [1] and verified in an exact nu- 
merical determination of ground states for interfaces in various 
types of disorders [2]. The main idea is to put the elastic system 
in a quadratic potential well, which acts as a large-scale cutoff 
and makes the problem well-defined. The shift between the 
center of mass and the center of the well is proportional to the 
renormalized force and its fluctuations are the quantity com- 
puted in the FRG [1]. The results of [2] show a remarquable 
agreement in the statics between the measured A (it) and the 1- 
and 2-loop predictions from the Functional RG [3-5]. These 
ideas and numerical tests have been extended to the depinning 
transition [6, 7] in the case of local elasticity, and to reaction 
diffusion models [8]. Finally, a first-principle calculation of the 
distribution of avalanches from the FRG was performed and 
verified by numerics [9-11]. An outstanding challenge is to 
test these predictions in experiments. 

The depinning of the contact line of a fluid on a disor- 
dered substrate has been studied experimentally [12-15] and 
its critical scaling established. Since gravity naturally creates 
a quadratic well, this raises the interesting possibility of mea- 
suring a FRG fixed-point function for depinning or, conversely, 
to learn more about the physical system using these (univer- 
sal) fluctuations as a new probe. To do so, the capillary length, 
which provides the well, does not need to be larger than the 
measurement scale. In fact, the finite capillary length is used as 
an advantage. 



Consider a fluid in a large reservoir and its contact line (CL) 
on a plate, parametrized by (x, u(x)) within the plate (fig. 1). 
Its energy can be modeled as 



n\u] 



dx 



[u(x) 



+ V(x,u(x))+£[u], (1) 



where £ [u] is a non-local elastic energy invariant under u(x) —> 
u(x) + const. The pinning force —d u V(x,u) is a local 
quenched random function. L is the total length of the con- 
tact line. Additional boundary terms, not written here, affect 




Fig. 1: Sketch of the experimental setup. The size of the image in the 
inset is 1.5 mm. 
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(degrees) 


(mN/m) 


(mm) (xl0~ 3 Pa.s) 


H 2 at 15K 


-40 


2.79 


1.94 0.021 


isobutanol 


-40 


20.9 


1.64 4.0 



Table 1 : Wetting properties 



the line near the boundaries x = and x = L. The position w 
of the center of the well with respect to the plate, which is also 
the equilibrium position of the CL without disorder, is fixed 
by the height in the reservoir far away from the contact line 
Woo, up to a constant shift w = + ^/2L C ^1 — sin 8, where 
Lc = V7/ ' 99 i s tne capillary length and 8 the contact angle. In 
the simplest model of the fluid surface the elastic energy £ [u] in 
Fourier space is £[u] = | J ^ekU-kUk , with an elastic kernel 
£fe = m 2 eg(k/kg). The curvature of the quadratic well is m 2 = 
jkg with kg = \pl sin(#) j(L c \J\ + sin 8). The scaling func- 
tion eg(x) is often approximated by 6^/2(2;) = \J x 2 + 1 — 1, 
but can be computed for any 8 [16]. 

One drives the system by slowly immersing the plate at a 
constant velocity v, hence w — vt. The center of mass 

u(t) ■= y [ dxu(x,t) (2) 
L Jo 

is fluctuating as a function of w, and contains valuable infor- 
mation about the system. In particular, the second cumulant 
A(w), defined as 

A(vt) := Lm 4 ([u(t + r) - v(t + t)][u{t) - vt]) c (3) 

with w — vt is exactly the renormalized disorder correlator 
computed by the FRG theory. Here and below (...) denote 
translational averages in the direction of the motion of the CL. 
Thus a direct measurement of A(w) is possible for the CL mov- 
ing on a disordered substrate. 

Two different systems have been used. The first one is liquid 
hydrogen on a cesium substrate, denoted H 2 /Cs hereafter. As in 
most previous experiments using HVCs, the Cs substrate is pre- 
pared at low temperature, yielding a rough surface with dense 
defects whose size is of the order of 10 nm [12, 13]. Such a 
small value leads to very small distortions of the CL and pre- 
cludes any optical observation. We have thus annealed the Cs 
layer up to 250 K. This creates a large-scale structure (possibly 
corrugation) which is visible but cannot be characterized since 
its typical length scale is still below the optical resolution. Its 
typical length scale is in the micrometer range. 

The second system is isobutanol on a silicon wafer, denoted 
iso/Si hereafter. We create a well-controlled disorder by pho- 
tolithographic technics: The Si wafer is decorated by random 
silanized square patches (10 x 10 /1m 2 ) [17] which cover about 
22% of the total area. As the silanized patches are less wettable 
than Si, each patch can pin the CL. 

The relevant parameters (advancing contact angle 8, liquid- 
vaper surface tension 7, capillary length Lc and dynamic vis- 
cosity i]) are listed in table 1 for both sytems. Note that contact 
angle and capillary length are similar. In both systems defects 
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Fig. 2: Height of the contact line h 2 L c i w ) averaged over 2L C , as 
a function of the position w of the plate (system: iso/Si). The fast 
depinning events (upwards) are clearly visible. Between them, the CL 
is adverted downwards at the plate velocity v (here 1 fim/s). The 
straight line is the reference level ho. 

are strong so that the Larkin length is set by the size £ of the 
defects; £ is macroscopic so that thermal activation is irrele- 
vant [12]. 

A sketch of the experimental setup is shown in fig. 1. The 
substrate is dipped into the liquid bath at constant velocity v. 
The CL is imaged with a standard progressive-scan CCD cam- 
era for velocities up to 10 fim/s. For the room-temperature 
system, we have also used a high-speed camera in order to anal- 
yse the CL dynamics up to v = 1200 fim/s. For each time t, 
the CL profile is digitized to obtain the CL height h(x, t). The 
correspondence with the above notations is 

h(x, t) — h = u(x, t) — vt — (u(x, t) — vt) . (4) 

To define the reference level h we first calculate the average 
height hi(t) over a CL length I: hi{t) := j L h(x, t) dx. An 
example is shown in fig. 2. Then h is obtained as the time- 
averaged value of hi(t). As the camera is fixed with respect 
to the liquid container, one expects h to be constant. How- 
ever, we allow for a slow drift of h for two reasons: First, as 
the volume of the container is finite, dipping the plate causes 
a slight increase of the asymptotic level which results in a lin- 
ear increase of the reference level ho- Secondly, we find that 
the slow variations of hi(t) are not linear, due to large-scale 
variations of the plate properties (defect size, mean contact an- 
gle, etc.). Hence, h (t) is determined experimentally by fitting 
the whole data set hi (t) by a polynomial function. As we are 
interested in correlations over distances two orders of magni- 
tude smaller than the total swept distance, this only shifts the 
correlator A by a constant. 

In the following, instead of h(x,t), we use h(x,w) := 
h(x,w/v) — h , where w = vt, making it easier to compare 
runs at different velocities v. We define the experimentally 
measured correlator A as 

A(w-w') := {hi(w)hi(w')) . (5) 



p-2 



Height fluctuations of a contact line: a direct measurement of the renormalized disorder correlator 



600 
500 



3. 



200 

<< 




/ ([Am) 




Fig. 3: Scaling of the correlator A(0) := (hi(w)hi(w)) as a function 

of the averaging length I along the CL (here iso/Si). The asymptotic Fig. 5: Inset: The disorder correlator A(w) for H 2 /Cs, with error 



scaling I 1 is marginally reached for I = L c 
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Fig. 4: Roughness W(l') of the contact line for iso/Si. Due to the lim- 
ited range between the defect size and the capillary length, no scaling 
is achieved. 



We need to choose the CL length I over which h is aver- 
aged. It must be larger than the capillary length L c , and as 
large as possible, while remaining notably smaller than the 
plate size L « 20 mm, since strong distortions occur at the 
edge of the plate. Moreover, slow changes of the plate prop- 
erties cause slow variations of the capillary rise on the plate. 
These are small compared to the fluctuations of the local CL 
height, but one expects that A varies like l^ 1 (for / ^> L c ). 
So, the larger I, the stronger the effect of the large-scale het- 
erogeneities. As a compromise, we have chosen I — 2Lq- 
As shown in fig. 3, A(0) has almost reached the asymptotic 
behavior for I — 2L C . Once I is chosen, an accurate deter- 
mination of A(w' — w) requires that the CL explore a large 
number of pinned configurations. In other words, the CL has 
to sweep a range of w much larger than the width of the line. 
Following [14, 15], we define the CL width W at a scale I' as 
W 2 (V) = ((h(x,t)~h(x + l / 1 t)) 2 ), where the average is taken 
over x and over all the successive configurations. For iso/Si, we 
find that the width W ~ 40 fim for V > 2L C ; for H 2 /Cs, the 
width is much smaller and difficult to measure, but we estimate 



bars estimated from the experiment. Main plot: The rescaled disorder 
correlator A(w) / A(0) (green/solid) with error bars (red). The dashed 
line is the 1-loop result from equation (6). 
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Fig. 6: Inset: The disorder correlator A(w) for iso/Si at v — 1/im/s 
up to w) = 35fj,m, and then at v — 10/im/s for w > 35/^m, with 
error-bars as estimated from the experiment. Main plot: The rescaled 
disorder correlator A(w)/A(0) (green/solid) with error bars (red). 
The dashed line is the 1-loop result from equation (6). 



that W ~ 5 /urn at large scale. The CL sweeps about 80 mm 
for iso/Si and about 5 mm for H2/CS. The plate displacement 
is thus much larger than W. This is necessary in order to get a 
reproducible shape for A. As a consequence, the disorder and 
the substrate properties have to be fairly homogeneous over an 
area of at least 1 cm 2 , which is quite difficult to achieve. In 
particular, we have first tried to use water on glass decorated 
with Cr defects as in [14], but we could not get rid of strong 
variations of the capillary rise. 

When comparing the two systems used in this work, the main 
difference lies in the defect size £, which is roughly one or- 
der of magnitude smaller for H2/CS compared to iso/Si. The 
drawback of the low-temperature system is the poor character- 
ization of the disorder which cannot be resolved optically. But 
the small value of £ is rather an advantage as the separation be- 
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tween the small-scale cut-off £ (defect size) and the large-scale 
cut-off L c (capillary length) is about 3 orders of magnitude for 
H2/CS, while only about 2 for iso/Si. As a consequence, it is 
difficult to observe the scaling regime for iso/Si. (For instance, 
W{V) does not exhibit scaling in fig. 4). The case of H2/CS 
is closer to previous experiments of the same type which did 
exhibit scaling. 

We now discuss our experimental results and their compari- 
son to theory. The raw data for the correlator A(w) defined in 
(5) are shown in the insets of figs. 5 and 6. To compare with 
theory it is useful to define a dimensionless disorder correla- 
tor as in [2] and plot the ratio Y := A(iu)/A(0) as a func- 
tion of the variable z with w = z J Q dw' A(u/)/A(0) so that 
the area under the curve is unity. The resulting function has 
no more free parameter and the functional RG theory of pinned 
systems predicts that in the limit of small m: (i) it should be the 
same function for all systems in a given universality class (UC); 
(ii) this function should exhibit a linear cusp near w = + . In 
figs. 5 and 6 we have plotted these dimensionless correlators 
(solid lines). Also shown in figs. 5 and 6 is the 1-loop predic- 
tion 

Z = V^i-ioop - 1 - In Yi -loop / J dy^y - 1 - Iny (6) 

for the universality class described by Eq. (1), i.e. quasi-static 
depinning with irrelevant non-linear terms, the only class for 
which this function has been computed analytically yet. One 
sees a very good agreement between data and this 1 -loop FRG 
prediction. A similar agreement was observed numerically both 
in the statics and the dynamics of pinned systems with local 
elasticity [2, 7]. At this stage we take this as a clear signature 
that we are dealing with a pinned system, and that a description 
using Eq. (1) is possible. A significantly higher precision and 
statistics (of a factor 1=3 10 — 20) would be required to reach the 
one which was achieved in numerics [2,7]. That would provide 
a decisive test on the universality class, and shed light on the 
debate about the large observed value for £ w 0.5 [15], while 
analytical predictions based on model (1) lie around £ = 0.4 
[4, 18]. Closer examination of the data in fig. 6 shows that devi- 
ations from Yi-ioop(z) can be mainly accounted for by rounding, 
which thus must be better controled and quantified. 

For iso/Si, we have performed the experiment at different 
velocities, in order to see the approach to the fixed point, and 
check for experimental artifacts. This is shown in fig. 7. We 
first comment on the shape of the correlator at the highest ve- 
locity (v = 1200 /im/s). The flattening at the origin is ex- 
pected as the system is driven away from the depinning thresh- 
old (i.e. away from the fixed point). However the large oscil- 
lations in A(w) (better visible in its derivative) are surprising. 
We believe them to be due to surface waves of the liquid, pos- 
sibly excited by the motor moving the plate. Decreasing the 
velocity, the cusp at the origin becomes more and more pro- 
nounced, though even at the lowest velocity v = 1 /im the first 
derivative A'(w) is not yet monotonically decreasing, as pre- 
dicted by theory for v = + . 

One possible origin of the observed rounding may be that the 
dissipation is not simply due to viscous shear in the meniscus 
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Fig. 7: The disorder correlator A(it) (top), and (minus) its first (bot- 
tom, main plot) and second (bottom inset) derivatives. The velocities 
are v — l^m/s (dashed-dotted/pink), v = 10/im/s (solid/black), 
v — 100/^m/s (blue/dashed) and v — 1200/im/s (green/dotted). Sta- 
tistical errors for A'(w) are ~ 5%. (The derivatives are obtained by 
fitting with a polynomial of degree 50 to 100 for < w < 120.) The 
non-zero value of A'(0 + ) is the signature of the predicted linear cusp. 



but involves complex microscopic processes at the solid sur- 
face. This is demonstrated by the fact that the behavior of the 
correlator for vanishing z is sensitive to the cleaning procedure 
of the plate, and by the fact that the homogeneous bare sili- 
con surface displays an intrinsic hysteresis (for a discussion, 
see [19]). At this stage, we have little idea how the underlying 
microscopic disorder could change the behavior of the CL at 
the scale of the macroscopically patterned defects, and how to 
predict the possible resulting rounding effect. 

On the other hand, there are also two fundamental reasons 
why the system could be slightly off the fixed point: (i) round- 
ing by the velocity, as observed above and expected form the 
theory (although most models predict that A'(0 + ) = at 
v > 0; this point is still debated [20]) (ii) as suggested by the 
width data in fig. 4, the ratio £/Lc ( or equivalently m) is yet 
too large to have reached the fixed point. The shape of A'(w) 
is qualitatively what is expected: A'(0 + ) is strictly positive, 
hence the system is above the Larkin scale, but it is signifi- 
cantly smaller than its putative fixed-point value, since A'(w) 
should be monotonically decaying there. Similarly A"(w) < 
near w — while the fixed-point value is expected to be pos- 
itive. This interpretation, if confirmed, is interesting, as it im- 
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Fig. 8: Dimensionless parametric plot of Q(w) / A(0) 2 versus 

1 — 4^ j for both experiments H2/CS (red/grey) and iso/Si at 

v — lOpim/s (black). The thin red (grey) line is the mean-field pre- 
diction. The thick dashed green (grey) line is an extrapolation based 
on a loop expansion, as discussed in the text. 
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Fig. 9: The measured dimensionless avalanche-size distribution func- 
tion p(s), as defined in the text, for iso/Si. For comparison an an- 
alytical estimate, defined in Eq. (12), based on a one loop calcula- 
tion with the scaled elastic kernel e n /2{x). The small-size cutoff is 
log 10 (s) = —2, corresponding to the size of a defect. 



plies that A (id) is a sensitive new probe, which can be made 
quantitative, to test how far the system is from criticality. 

More information can be obtained from the experimental 
data by computing the third cumulant: 



S 3 (w-w') ■= ([hzH-hzK)] 1 



It is convenient to plot the function 
1 



Q(w) := I J dw'S 3 (w') , 



Indeed, the FRG predicts that 

Q(w) 



A(0)* 



A 



A(w) 



(7) 



(8) 



(9) 



with an exact equality and the universal amplitude A — 1 at the 
mean-field level (i.e. resummation of tree diagrams, d > d uc = 
2). The proportionality (with no attempt to measure A) was 
checked numerically [7] for depinning with local elasticity. In 
fig. 8 we plot both sides of eq. (9): we see that the proportion- 
ality holds very nicely and that the two experiments fall on top 
of each other, confirming the universality of the slope A, which 
is measured to be A ps 0.53. The FRG calculation, using the 
scaled elastic kernel 6^/2(2;) yields A = 1/(1 + |e) + 0(e 2 ) 
with e = 2-d (see (E12) of [10]) which yields 9/17 = 0.53 
and 0.11 for the two Pade appro ximants at e = 1, while the 
kernel eg(x) with 6 = 40° gives [16] A = 0.48 ± 0.13 (shown 
in fig. 8). These values are reasonable, given that the correction 
to mean field is large and one cannot hope for high precision. 
Also note that deviations from the functional form in Eq. (9) 
are expected at 1-loop order, but they should be small as found 
numerically in [7]. We conclude that the agreement between 
experiment and theory is satisfactory for the third cumulant. 

The properties of the CL can also be characterized by the 
distribution of the sizes S of avalanches, or forward jumps 
(see fig. 2), where by definition S is the area swept during the 



avalanche. At the critical point, m = 0, the distribution is ex- 
pected to be a power law characterized by an exponent r. At 
small m > 0, the correlation length is finite and the avalanche- 
size distribution P(S) for S ^S> S m \ n is cut off at scale 



Q — 



2(S) 



(10) 



where here and below (. . .) denotes an average over P(S). One 
expects that in the variable s := S/S m the avalanche-size dis- 
tribution exhibits universality, i.e. independence of short scales, 
and that, for 1 < r < 2, it takes the form 



P(S)dS:=£>p(£)^ 

\ &m / ^ in 



(ID 



The function p(s) is universal, and by construction from (10) 
and (11) normalized s.t. L ds sp(s) = 1 and L ds s 2 p(s) = 
2 [21]. The analytical prediction [10] forp(s) reads, based on 
the model (1) and with the scaled elastic kernel e n / 2 (x): 



p(s) = A's T exp 



B' 



D'i 



(12) 



with B' = 1 + A( 7 - 2)e, C = ^p, D' = A' given by 
the normalization and the exponent 8' = 1 + 1. This prediction 
is exact to first order in e = 2 — d, and to produce our analytical 
estimate we set e = 1 and rescale both axis to ensure the two 
normalization conditions. For the r exponent, we have used 
the conjectured relation t(£) = 2 — l/(d + £) since we proved 
[10] that the latter is exact at least to one loop. Inserting the 
measured («0.5 yields the prediction r w 4/3. 

The size distribution is measured only for iso/Si with a stan- 
dard camera (acquisition rate: 15 Hz), at a velocity v = 
1 /jm/s. We consider that an event occurs if the local displace- 
ment between two successive images is larger than a threshold 
6h ~ 8 /jm. This threshold is about the defect size, which is the 
smallest displacement occuring in an avalanche, and is about 10 
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Fig. 10: Z(A) defined in (13), from bottom to top: the mean-field 
prediction (red, solid), the experimental data (black dots), and the ex- 
trapolation based on the one loop calculation for a scaled elastic kernel 
eg(x) with 6 = 40° (green, dashed). 



moment of the avalanche-size distribution can similarly be re- 
lated to the third cumulant of the center-of-mass fluctuations: 
(S 3 ) (S) /(3{S 2 ) 2 ) = A, with A denned in eq. (9). From 
p(s) we find A = 0.77, while the relation between Q(w) and 
A(w) yields A = 0.53. The agreement in both cases is only 
fair and more experiments are needed. In particular, a smaller 
ratio £ /Lq would be helpful. 

To conclude: By examining the fluctuations of the mean 
height of the contact-line at depinning, we have measured the 
renormalized disorder-correlator A (w). The latter is the central 
object of the functional RG field theory, and its predicted cusp, 
which is the sign of metastability, shocks and avalanches, was 
under intense debate from the field theory side. Here we have 
made the first comparison between experiment and theory. It 
shows qualitatively and quantitatively that the ideas in the lat- 
ter are correct, and opens new ways of quantifying the former, 
calling for new experiments. 

We thank G. Borot and M. Pettersen for help with the exper- 
iments and A. Rosso for discussions. 



times larger than the resolution on the CL position. The area S 
swept in a single avalanche is then computed as the sum of the 
displacements occuring before the CL is pinned again. We de- 
tected 16 x 10 3 avalanches for a swept area of 3.5 x 20 mm 2 . 
We have checked that the resulting distribution is not sensitive 
to the value of the threshold and that the acquisition rate is fast 
enough to avoid lumping together uncorrelated avalanches. 

Experimental data and predictions are shown in fig. 9. We 
see that the parameter-free scaling collapse is quite good down 
to s = 10~ 2 . This value of s corresponds to S = 100 /mi 2 , 
which is the cutoff expected from the size of the defects whose 
area is 100 /mi 2 . However our accuracy is limited by the 
available sizes, and e.g. not sufficient to discriminate between 
r(C = 0.5) = 4/3, and r(C = 0.4) = 1.29. While p(s) has 
strong statistical fluctuations, the characteristic function 



Z{\) 



dsp(s) 



1 



(13) 



can be measured quite accurately as shown in fig. 10. The 1- 
loop extrapolation (e = 1) is closer to the data than the mean- 
field result and using the elasticity €9=40° (x) (shown in fig. 10) 
is also closer than using ie=n/2{x) (not shown). Hence the 
comparison to field-theory predictions [10, 16] is satisfactory. 

We now come to a final test of the avalanche picture under- 
lying the FRG calculations. According to [9, 10], the slope 
A'(0 + ) of the linear cusp is proportional to the scale S m of 
the avalanche-size distribution, 



L 



A'(0+) 



_ 9 = <^!> 

~ m ~ 2 (S) 



(14) 



We find S m ~ 9000/mi 2 . At the lowest velocity, |A'(0+)| = 
1.75/mi 2 andi = 3500/im, one finds L|A'(0+)| ~ 6100 /im 2 . 
However, estimating A'(0 + ) is difficult because we do not un- 
derstand the origin for the rounding of A(w) when w — ► 0. 
An upper bound for |A'(0 + )| is the slope of A' at the inflec- 
tion point, which yields L|A'(0+)| < 14000 /mi 2 . The third 
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